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Abstract. The metric dimension of a graph G is the size of a smallest 

subset L C V(G) such that for any x,y £ V{G) there is a z £ L such 

that the graph distance between x and z differs from the graph distance 

fSJ between y and z. Even though this notion has been part of the literature 

T-H for almost 40 years, the computational complexity of determining the 

^^ metric dimension of a graph is still very unclear. Essentially, we only 

vN know the problem to be NP-hard for general graphs, to be polynomial- 

r^ time solvable on trees, and to have a log n-approximation algorithm for 

D general graphs. In this paper, we show tight complexity boundaries for 

pLn the Metric Dimension problem. We achieve this by giving two com- 

(-^ plementary results. First, we show that the Metric Dimension prob- 

^sj lem on bounded-degree planar graphs is NP-complete. Then, we give 

a polynomial-time algorithm for determining the metric dimension of 
outerplanar graphs. 



1 Introduction 



Given a graph G — iV^E), its metric dimension is the cardinality of a smallest 

set L C y such that for every pair x^y €V , there is a z G L such that the length 

vj of a shortest path from z to a; is different from the length of a shortest path from 

vQ z to y. In this case we say that x, y are resolved by z and L. Elements of the set 

iy~j L are called landmarks. A set L C V that resolves all pairs of vertices is called 

CNJ a resolving set. The problem of finding the metric dimension of a given graph G 

Cn is called Metric Dimension, but is also known as Harary's problem, and the 

[~^ locating number or rigidity problem. The problem was defined independently 

^D by Harary and Melter [T^ and Slater [TH], and has received a lot of attention 

^~~' from researchers in different disciplines (see e.g. |2|13j and references therein). 

. . For example, a recent survey by Bailey and Cameron [2] notes an interesting 

^ connection to group theory and graph isomorphism. 

Computing the metric dimension has proved to be a hard problem, and there 
are very few results on the computational complexity of the problem. The de- 
ed cision version of Metric Dimension is known to be NP-complete on general 

graphs [11|15I16] . KhuUer et al. [T6j gave a linear-time algorithm to compute 
the metric dimension on a tree. Characterizations for graphs with metric di- 
mension 1 and 2 are also given in the same work. Moreover, |16, also provides 
a 2 log n-approximation algorithm for the metric dimension of any graph. Beer- 
liova et al. [3 proved that the metric dimension of a general graph cannot be 
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approximated within o(lnn), unless P=NP. Hauptniann et al. [13j improved on 
this by giving a (1 — e) logn inapproximabiUty result under mild complexity as- 
sumptions. They also show that determining the metric dimension is APX-hard 
for bounded-degree graphs, and APX-complete for superdense graphs (graphs 
for which the complement is a sparse graph). To the best of our knowledge, no 
further results are known about the complexity of this problem. 

Our Results In the present work, we narrow the tractability gap of Metric 
Dimension. From the hardness side, we show that Metric Dimension on pla- 
nar graphs, called Planar Metric Dimension, is NP-hard, even for bounded 
degree planar graphs. Even though Garey and Johnson considered Metric Di- 
mension over 30 years ago in their book [11, Problem GT61], but the complexity 
of the planar version of the problem remained unknown thus far. From the al- 
gorithmic side, we show that there is a polynomial-time algorithm to find the 
metric dimension of outerplanar graphs. This significantly improves the 15-year- 
old algorithm on trees by KhuUer et al. |TB]. 

Our algorithm for outerplanar graphs presents a leap forward in the type of 
problems that can be solved on such graphs. Many graph problems, e.g. inde- 
pendent set, are local in the sense that a vertex has direct effect on its neighbors 
only — indirectly the effect propagates farther. The main obstacle in the Met- 
ric Dimension problem is that it is nonlocal: a landmark can resolve faraway 
vertices. This means that standard approaches fail to work for it. In her seminal 
paper, Baker writes that the dynamic programming technique therein "wi// work 
for problems that involve local conditions on nodes and edges" [41. This does not 
include Metric Dimension. For the same reason, the generalizations of Baker's 
work by Eppstein [10 and Courcelle [6 do not apply either, and a different kind 
of approach is necessary. 

The crux then to both our hardness and algorithmic results is our ability 
to constrain the effects of a landmark to a small area. The NP-hardness proof 
constructs a specific family of planar graphs for which Metric Dimension is 
essentially a local problem. The algorithm uses a tree structure to traverse the 
graph, together with several data structures that track the influence of landmarks 
on other vertices. The structures that we employ differ significantly from the ones 
used in standard bounded-treewidth techniques. Moreover, whereas outerplanar 
graphs have constant treewidth [5] , the tree structure used in our approach leads 
to a decomposition that can have arbitrary width. 

We believe that our algorithmic techniques present a generalization of Baker's 
work, and are of independent interest. In particular, they might lead to (new) 
algorithms for a broad class of non-local problems (see also Section [4]). 

Overview of the NP-Hardness Proof As a corollary of the work by Dahlhaus et 
al. i8], we prove a new version of Planar 3-SAT to be NP-complete. We reduce 
this problem to Metric Dimension. This is done by constructing a planar 
graph consisting of clause gadgets and variable gadgets. Let n be the number 
of variables. Each variable gadget must have 4 landmarks, 3 at known, specific 
locations, but for the fourth we have three different choices. They correspond 
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to the variable being true, false, or undefined. These An landmarks are a metric 
base if and only if they resolve all pairs in clause gadgets, which happens only 
if they correspond to a satisfying truth assignment of the SAT-instance. 

Overview of the Algorithm First, we characterize resolving sets in outerplanar 
graphs by giving two necessary and sufhcient requirements for an arbitrary vertex 
set to be a resolving set. Then, taking as a base the dual of the biconnected 
components of the graph G, we define a tree T. Vertices of T corresponds to 
faces and cut vertices of G, and edges to inner edges and bridges of G. The 
algorithm uses dynamic programming to process T, starting at the leaves and 
advancing towards the root. Data structures called boundary conditions and 
configurations are used on edges and vertices of T, respectively, to describe how 
the landmarks are positioned. 

Boundary conditions track the effects of landmarks placed in the already 
processed part of the graph and the possible effects of sets of landmarks to 
placed in the unexplored paths of the graphs. 

The main algorithmic novelty lies in the configurations, which control the 
process of combining the boundary conditions on edges towards children of the 
current vertex v' € V{T) into a boundary condition on the edge towards the 
parent of v' . The configurations highly depend on the vertices of G represented 
by v' . However, even though the number of vertices of G represented by v' might 
be unbounded, we can show that the total number of relevant configurations is 
only polynomial. The use of configurations also presents a stark contrast with 
the techniques used in bounded treewidth algorithms, where the combination 
process commonly is a simple static procedure. 

Preliminaries For basic notions and results in graph theory, we refer the reader 
to any textbook on the topic, e.g. [5j- AH graphs are finite, undirected, and unless 
otherwise stated, connected. We use the notation (u, v) to denote an edge from 
u to V. By distance d{u,v) we mean the graph distance. The vertex and edge 
sets of H are denoted by V{G) and E{G), respectively. Recall that a graph is 
planar if and only if it does not contain a subgraph that is a subdivision of K^ or 
K33, and it is outerplanar if and only if it does not contain a subdivision of K4 
or K2^3. A graph G has a cut vertex if the removal of that vertex disconnects the 
graph into two components. A graph is a biconnected if it has no cut vertices. If 
G is a biconnected outerplanar graph, G has a planar embedding in which the 
edges on the boundary of the outer face form a Hamiltonian cycle. We call those 
edges outer edges and the rest of the edges are called inner edges. Given G and 
V G V{G), Af{v) denotes the set of neighbors of v in G. Given a set S, we denote 
by V{S) the power set of S. 

2 NP-hardness on planar graphs 

It is well known that 3-SAT is NP-complete (see [H]). We require a special 
planar version of 3-SAT to be NP-complete. 
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Definition 1 Let ^ be a boolean formula, which uses the set of variables V and 
has the set of clauses C. Then the graph G^ = {V yj C, E), where E — {(v, C) \ 
V €V,c€ C,v ^ C}, is the clause- variable graph of^. 

With the notation w e C we mean that variable v (or its negation) occurs in 
clause C. Observe that Gq, is always bipartite. 

Theorem 2 ([8, p. 877]) The problem of deciding whether a boolean formula 
\P is satisfiable is NP- complete, even if 

— every variable occurs in exactly three clauses (twice positive, once negative), 

— every clause contains two or three distinct variables, and 

— Gip is planar. 

As a corollary of Theorem [2] we get the following result, which is the starting 
point of our work. 

Corollary 3 The problem of deciding whether a boolean formula W is satisfiable 
is NP-complete, even if 

— every variable occurs exactly once negatively and once or twice positively, 

— every clause contains two or three distinct variables, 

— every clause with three distinct variables contains at least one negative literal, 
and 

— G^ is planar. 

We call this decision problem 1-Negative Planar 3-SAT. 

Proof. Let tf' be a boolean formula satisfying the constraints of Theorem [2j By 
modifying tf' we will construct a formula 'F' that fulfills all the constraints of the 
theorem statement and is satisfiable if and only if if' is satisfiable. 

We only need to eliminate those clauses containing three positive literals. 
Suppose that a; V y V z is such a clause of ^ with distinct variables x, y, z. Now 
add a new variable x' , and replace the original clause by the clauses xV x' and 
-la;' y yy z. This completes the construction of If"'. As this construction replaces 
some edges of G^ with paths, it preserves planarity. 

Given a satisfying truth assignment of W, we get a satisfying assignment 
of S'' by setting x' = —■x. A satisfying assignment of W' implies a satisfying 
assignment of \I/. So \lf' is satisfiable if and only if W is. The theorem now follows 
straightforwardly from Theorem [2] 

To prove that Planar Metric Dimension is NP-hard, we wiU give a reduction 
from 1-Negative Planar 3-SAT. The idea behind the graph constructed in 
this reduction is the following. Given an instance ^ of 1-Negative Planar 3- 
SAT, we first find a planar embedding of its clause- variable graph G^f. We then 
replace each variable vertex of G^ by a variable gadget (see Figure [l]) , and each 
clause vertex of G^p by a clause gadget (see Figure^. By identifying vertices 
of variable gadgets and vertices of clause gadgets in an appropriate way (see 
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Fig. 1. The variable gadget. 
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Fig. 2. The clause gadget for a clause with two variables (left), and the gadget for a 
clause with three variables (right). 
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Figure I4]), we obtain a planar graph Hq^ that wiU be our instance of Planar 
Metric Dimension. 

We now describe our construction in detail. Consider a planar embedding 
of G,p, which can be found in linear time 14 . We first replace each variable 
vertex of Ci: by a variable gadget. In Figure [l] the vertices marked by unfilled 
circles will be identified with vertices from a clause gadget later on. There are 
three groups (connected components) of such vertices in the figure. The groups 
containing vertices (ii,/i) and (is, /a) will be identified with vertices in clause 
gadgets where this variable appears positively in the corresponding clause; the 
group containing (^2, /2) will be identified with vertices in clause gadgets where 
this variable appears negatively. By rotating and contorting the variable gadget 
appropriately, we can ensure that the three groups point into the right direction 
(i.e. the negative-appearance group faces the clause vertex where the variable 
appears negatively). 

Next, we replace the clause vertices by clause gadgets. The exact gadget we 
use depends on whether the clause contains two or three variables (see Figure [2]). 
We restrict our description to the three- variable case, as the two- variable case is 
similar and simpler. In Figure [2J the vertices marked by unfilled circles will be 
identified with vertices from a variable gadget. There are again three groups of 
such vertices, one for each variable occurring in the clause. 

Obviously, we will identify the i-vertex of a variable group with the i-vertex 
of a clause group, and the same for the /-vertices. We call this matching. 
It is not entirely straightforward to do this matching in a manner that pre- 
serves planarity. Consider the way in which the groups and the t and / ver- 
tices appear on the boundary of the clause gadget. In Figure [2] the pairs ap- 
pear in order (t, /), (t. /), (/, i) clockwise starting from the top. As illustrated 
in Figure [SJ (t, /), (/, t), (/, i) is also possible. The remaining two alternatives, 
(^i /); (^7 /)i (^i /) and (/, t), (/, t), (/, t) are to be avoided. This is accomplished 
by choosing a variable appearing negatively in the clause and mirroring the cor- 
responding variable gadget around the axis Ti — F (see Figure IT]). This does not 
affect our ability to connect the variable to other clauses. 





h h 



Fig. 3. Alternative planar embeddings of the the clause gadgets. 



On the Complexity of Metric Dimension 7 

This completes the construction. Call the resulting graph H^p, which is planar 
and has bounded degree by construction. We remark that each variable appears 
once negatively in ^, and once or twice positively. So if the variable appears only 
twice, then (ii,/i) or it^^f^) in the corresponding variable gadget will not be 
identified with a group of vertices in a clause gadget. 

In Figure [4] we can see an example of the reduction and the resulting planar 
graph from the specific instance of 1-Negative Planar 3-SAT (^xi V a;2) A 
(a;i V ^X2 V x-i) A (x2 V 0:4 V ^2:3) A (^^4 V x^). 

We now make several observations about the graph Hxp and the way land- 
marks need to be positioned on it. 

Each /-vertex is contained in a triangle, say with other vertices r, s. Observe 
that r and s can only be resolved by a landmark on r or s. We call these forced 
landmarks. In fact, in any smallest set of landmarks exactly one of r, s will be 
a landmark. Then it follows by construction that H4, requires exactly 3n forced 
landmarks, where n is the number of variables of W. 






Fig. 4. The planar graph obtained for (^a;ia;2)(a;i^a;2a::3)(a;2a::4^X3)(^X4a;3) 



Using the forced landmarks, we can resolve most pairs of vertices, as shown 
by the following lemma. We say that Ti, T2, Ni, N2 and F are strictly inside 
the variable gadget. 
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Lemma 4 Letx,y be a pair of vertices, which is not {wi,W2}, {Ti,T2}, {Ti, A^i} 
or {T2, Ni} from the same gadget. Then the pair is resolved by forced landmarks. 

Proof. It is straightforward to check the cases where both x and y are in the 
same clause or variable gadget. There are two remaining cases: either x, y are in 
different clause gadgets, or x is strictly inside a variable gadget and y outside 
that gadget. 

Consider the first case, that is, x,y are in different clause gadgets. Denote 
the gadget containing x hy g^, and the gadget containing y by gy. Let z^ be 
a forced landmark that is closest to x, and let Zy be a forced landmark that is 
closest to y. Without loss of generahty, d{x,Zx) < d{y,Zy). We will show that 
d{zx,x) < d{zx,y). Since x,y are in distinct clause gadgets, the shortest path 
P from y to z^ crosses at least one variable gadget g, and enters it through a 
group (i, /). Let w be the first vertex of {t, /} that occurs on P, and let zj be 
the forced landmark in the triangle connected to f.liw = t, then P contains at 
least two edges in g, and ii w = f, P contains at least one edge in g. In either 
case P also contains at least one edge in g^. Since d{zf,f) = 1, d{zf,t) = 2, in 
both cases the inequality d{w,Zf) < d{w,Zx) holds. Hence d{x,Zx) < d{y,Zy) < 
d{y, Zf) < d{y, w) + d{w, zj) < d{y, w) + d{w, z^) = d{y, z^). 

Now consider the second case, and assume that x is strictly inside the variable 
gadget in Figure [5] If y is in the picture, it can be readily verified that x and 
y are resolved by the forced landmarks of the variable gadget. We claim that 
if y is outside of the picture, then d{zi,y) + d{z3,y) > 7, where zi and Z3 are 
the forced landmarks in the triangles attached to /i and /s respectively. This 
implies that zi or Z3 is at distance at least four from y, whereas the distance of 
zi and Z3 to x is at most three, implying that x and y are resolved. 

To prove the claim, note that if shortest paths from y to zi and Z3 both 
contain /i (or equivalently /s), then d{zi,y) + d{z3, y) = d(zi, /i) + d(z3, /i) + 
2d(/i, y) = 4 + 2d{fi,y) > 6. Now consider the case where shortest paths from y 
to Zi , Z3 include /i , f^ respectively. Any path from /i to f^ that is not contained 
in the figure must cross at least two clause gadgets and one variable gadget, so it 
has length at least 5. This gives d{y, fi)+d[y,fz) > 5, and d{zi,y)+d{z3,y) > 7. 

It remains to analyze how the pairs excluded in Lemma |4] can be resolved. This 
will rely on the satisfiability of 'F of course (as described below) , but the following 
auxiliary lemma is crucial. 

Lemma 5 All pairs of vertices that are strictly inside a variable gadget are 
resolved if and only if there is a landmark strictly inside the variable gadget. 

Proof. It is easy to check that a landmark that is strictly inside a variable gadget 
together with the forced landmarks resolves all pairs of vertices that are strictly 
inside the gadget. If no landmark is strictly inside the variable gadget, then 
a shortest path from any landmark z to Ti or T2 contains ti or ^3. But then 
d{z,Ti) = d(z,T2). 

This lemma and the forced landmarks together imply that Hq, has metric di- 
mension at least 4n. With this fact in mind, we present the proof of our main 
result. 
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Fig. 5. The subgraph in Lemma |4] 



Theorem 6 Planar Metric Dimension is NP-complete, even on bounded- 
degree graphs. 



Proof. Let ^ be an instance of 1-Negative Planar 3-SAT with n variables. 
Construct the graph H^ in the manner described before. Constructing H,i, clearly 
takes time polynomial in the number of variables and clauses of ^Z^. We now claim 
that H^ has metric dimension An if and only if !f' is satisfiable. 

Assume that a satisfying truth assignment for •Z' is given. Place 3n forced 
landmarks. If a variable has value true, place a landmark on Ti in the corre- 
sponding gadget; otherwise, place a landmark on F. After applying Lemmata|4] 
and [5] we only need to check that pairs wi,W2 are resolved. But each such 
pair is resolved by the landmark strictly inside the variable that satisfies the 
corresponding clause, and we are done. 

Now assume that the metric dimension is in. We will construct a satisfying 
assignment for ^. Each variable gadget contains exactly one landmark, which is 
on Ti, Ni, or F. If the landmark is on Ti, set the variable to true. If the landmark 
is on F, set it to false. Otherwise the variable can be arbitrarily set to either 
true or false. It remains to show that because the pairs wi , W2 are resolved, 
the truth assignment is satisfying. Note that a landmark z resolves pair wi , W2 
if a shortest path from a landmark to either of them enters the clause gadget 
through some ti. If a shortest path from landmark z to wi or 1x12 intersects more 
than one clause gadget, it leaves the first clause through an /-vertex, after which 
it enters all subsequent ones through an /-vertex. But then wi , 1x12 in the final 
clause are not resolved. It follows that a landmark z resolves wi and W2 only if 
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it is in an adjacent variable gadget and the corresponding variable satisfies the 
corresponding clause. 

Following the claim, the reduction should construct H^ as the graph for the 
instance of Planar Metric Dimension and set k to An. 



3 Algorithm for Outerplanar Graphs 

In this section, we prove that Metric Dimension can be solved in polynomial 
time on outerplanar graphs. So let G be an outerplanar graph, given together 
with an outerplanar embedding. Note that a metric base of a disconnected graph 
is the union of metric bases of its component (with one exception: isolated ver- 
tices, an edgeless graph of n vertices has metric dimension n — 1). So we can 
safely assume that G is a connected graph. We can also assume that it has at 
least three vertices. 

As a technical trick we sometimes treat the midpoint of an inner edge e = 
ivi,V2) € E{G) as an actual vertex. The distances from this midpoint ver- 
tex Ve are such that d{ve,vi) = d{ve,V2) = \ ^'^d dlve^x) = min((i(we, wi) + 
d{vi,x),d{ve,V2) + d{v2,x)). 

3.1 Characterization of Resolving Sets in Outerplanar Graphs 

We first give a characterization of the effects of resolving sets in vertices and faces 
of outerplanar graphs. To this end, we introduce some notation. A bifurcation 
point associated with z, x, j/ is a vertex v farthest from z such that v is on shortest 
path from z to both x and y. More formally, w is a bifurcation point if it is on 
shortest paths z ^^ x, z ^^ y, and if any shortest paths v -^ x, v ^^ y intersect 
only in v. Notice that in an outerplanar graph the bifurcation point for each triple 
of vertices is unique. We define the function g : V{G) x V{V{G)) -^ V{V{G)) as 

g{v, L) = {w e Af{v) : d{z, w) = d{z, v) + 1 for aU z e L}. 

In other words, a neighbor w of w is in g[v,L) if for every z € L, v is on some 
shortest path z ^^ w (but not necessarily on every shortest path.) 

Any pair x,y € g{v, L) is left unresolved by landmarks in L. So any resolving 
set L satisfies the following: 

Requirement 1 Any vertex v G V{G) must have \g{v,L)\ < 1. 

If G is a tree, then the requirement is sufficient. The following lemma gives an 
alternative, simpler correctness proof for the algorithm by Khuller et al. |16j . 

Lemma 7 If G is a tree with at least three vertices, then a set L C V{G) is a 
resolving set if and only if it satisfies RequirementY^ 

Proof. We have already seen that any resolving set L satisfies Requirement [TJ 
Now assume that Requirement [l] is satisfied. We pick any two vertices x,y G 
V{G) and show that they are resolved. 



On the Complexity of Metric Dimension 11 

Since G has at least 3 vertices, there is at least one vertex v G V{G) with 
degree at least 2. Since \g{v,L)\ < 1 < |A/'(u)|, L is not empty. 

Choose any landmark z € L. li z resolves x,y, we are done. Otherwise, let v 
be the bifurcation point associated with z,x,y, and let Ui,'y2 be the successors 
of V on the shortest paths v -^ x,v -^ y (see Fig. p|. Since d{z,x) — d{z,y), 
we have d(y,x) = d{v,y). By assumption, g{v,L) can not contain both vi and 
V2- Without loss of generality vi ^ g{v,L). Then there is a landmark Z2 with 
d{z2,vi) < d{z2,v), and furthermore d{z2,x) < d{z2,y). □ 



Fig. 6. Proof of LemmalT] 



As stated earlier, the major difficulty of the metric dimension problem is that 
it is non-local. This is why Lemma [7] is useful. Although stopping short of giving 
an actual local characterization of resolving sets, it does make the effects of a 
resolving set local enough to devise a polynomial-time algorithm for trees. 

Our algorithm relies on a generalization of LemmalTlto outerplanar graphs. In 
this case Requirement [l] no longer implies that L is a resolving set. For example, 
if G is an even cycle and L contains two antipodal vertices, then Requirement [T] 
is satisfied, but L is not a resolving set. Therefore we need another requirement 
that provides an additional condition for the effects of a resolving set on the faces 
of outerplanar graphs. First, we need some auxiliary definitions and lemmas. 

Definition 8 Let z e V (G) he a landmark, and let C he either a single edge or a 
cycle — in particular, it may he a face. The representative of z on C is the element 
of V(C) closest to z, if it is unique. Otherwise outerplanarity implies that there 
are two closest vertices, which are adjacent. In this case the representative is the 
midpoint of those two vertices. 

We can make the following observations. Let C, C" be cycles such that 
V{C') C V{C) and C" is a face. If two vertices have the same representative 
on C, then they must have the same representative on C" as well. 

Lemma 9 Let G he a graph, and let x,y,z,Z2 € V{G). If neither z nor Z2 
resolves the pair x, y and there exist intersecting shortest paths z ^^ x, Z2 ^^ y, 
then a bifurcation point of z,x,y is also a bifurcation point of Z2,x,y. 
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Fig. 7. Requirement [2] 

Proof. Assume that w is a vertex which is both on some shortest path z '-^ x 
and some shortest path Z2 ^^ y- Let w be a bifurcation point of w,x,y. Then 
there exist shortest paths z-^ x, Z2'^ y both containing both v and w. We need 
to show that v is also in the intersection of some pair of shortest paths z "^ y, 

Zi "^ X. 

The proof rehes on the fact that the equaUty d{a, c) = d{a, h) + d{b, c) holds 
if and only if 6 is on a shortest path from a to c. 

By the triangle inequality and since neither z nor Z2 resolve a;,?/, 

d{z, y) + d{z2, x) < {d{z, v) + d{v, y)) + {d{z2.,v) + d{v, x)) 

— {d{z, v) + d{v, x)) + {d{z2,v) + d{v, y)) 
^d{z,x) + d{z2,y) 

— d{z,y) + d{z2,x). 

Therefore equality holds in the triangle inequality, and the claim follows. D 

We are now ready to present the other necessary requirement. 

Requirement 2 Let v' be a face of an outerplanar graph G and let L C V{G) 
have exactly two representatives zi,Z2 on v' . Let z^ and z^ he the landmarks 
with representative zi which occur first and last on the walk along the outer face 
starting at Z2, and define Z2, ^2 analogously. Assume that neither x^ nor X2 
resolves x, y, and that shortest paths z{ -^ x, z^ '^ y do not intersect. Let v be 
the bifurcation point of zf ,x,y, let u be the bifurcation point of x,zi,Z2, etc., 
as in Fig. [^ By assumption, v ^ t. Therefore the shortest paths s '^ t, t '^ u, 
w^ V, V '^ s form a cycle C (see again Fig. MJ. If v' is inside the cycle C , then 
one of z[, z\ must resolve the pair x,y. 

Note that the representatives of z( and Zj on C are v and t, respectively. 

The ass umption that G is outerplanar is essential for Requirement [2] and 
Lemma 10 below. In particular, the definition of z(, Z2, Z2, z\, as well as the use 
of representatives in the proof of Lemma |10[ relies on outerplanarity. 



Lemma 10 // G is an outerplanar graph, then any resolving set L C V{G) 
satisfies Requirement^ 
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Fig. 8. Neither of 21, z[ has representative v on C. 

Proof. Assume that L is a resolving set, and that v' , x, y are as in Requirement [2J 
Without loss of generality, there is a ^i S L with representative zi on v' such 
that d{zi,x) < d{zi,y). Now define z(,Z2,z{ as in Requirement pi Assume that 
neither z( nor 23 resolves x and y. Then we can define v, t, C as in Requirementbl 
Observe that the vertices v and t are the only vertices of C that do not resolve 
X and y. Therefore, while landmark zj has representative v on C, zi does not. 
Because of the location of z(, zi, z\ on the outer cycle, the representative of z\ 
on C is not v either (see Fig. ^ . The representative of z\ on C cannot equal t, 
as z\ and Zj have distinct representatives on v' . Therefore the representative of 
z\ resolves x and y, and so does z\ itself. D 

Now we are ready to generalize Lemma [7] to outerplanar graphs. This is 
a crucial result, since it characterizes resolving sets in a manner that allows 
dynamic programming. 

Theorem 11 If G is an outerplanar graph, then a set L C V{G) is a resolving 
set if and only if it satisfies Requirements^^ and[^ 

Proof. We have already seen that any resolving set satisfies the requirements. So 
assume that L C V{G) satisfies the Requirements, and choose any x,y £ V{G). 
We show that there exists a z G L that resolves the pair x, y. 

As in Lemma [71 L is non-empty. Choose z G L arbitrarily. If z resolves x 
and y, we are done; so assume that it does not. As in Lemma [7J let v be the 
bifurcation point of z,x,y, and let ui,i'2 be successors of v on some shortest 
paths V ^^ x,v '^ y respectively. By Requirement [T] there is a. Z2 & L such that, 
without loss of generality, d{z2, vi) < d{z2, v). If Z2 resolves x and y, we are done. 
So assume otherwise. 

If there exist intersecting shortest paths z '-^ x, Z2 ^^ y, then by Lemma [9] 
i; is on a shortest path 22 ^^ x. Then vi is also on such a path, and d(z2,vi) = 
d{z2 ,v) + l, a contradiction. Therefore no such pair of intersecting shortest paths 
exists. Define v, t, C as in Requirement [2] and let v' be a face inside C. If there 
exists a Z3 e L whose representative on v' is distinct from the representatives of 
2i and Z2, then its representative on C is neither v nor t. Hence 23 resolves x, y. 
If such a landmark 23 does not exist, then Requirement [2] implies the claim. D 



14 



Josep Diaz, OUi Pottonen, Maria Serna, and Erik Jan van Leeuwen 



3.2 Generalized dual tree 

We now introduce first of the the data structures that we need in our algorithms. 
The weak dual of an outerplanar graph G is a graph which has the faces of G, 
except the unbounded outer face, as vertices. Two faces are adjacent if they 
share an edge. The weak dual of a biconnected outerplanar graph is a tree. 



v„ ♦ 





C(Fi,i.i) 



B(i.i,F,) 



Fig. 9. An outerplanar graph G and its generalized dual tree T. The blue dashed line 
indicates how G is divided to B{Fi,vi) and B{vi,Fi). The sets on the edges of T 
indicate the values of s(e'). 



The generalized dual tree T = {V' , E') of an outerplanar graph G is defined 
as follows. T contains the weak dual of G, and also contains the subgraph of 
G induced by all cut vertices and vertices of degree one. For any cut vertex 
contained in a biconnected component, there is an edge from the vertex to an 
arbitrary incident face of the component. Observe that the resulting graph is a 
tree (see Fig. [o] for an example) . According to this definition, a cut vertex is a 
vertex of both G and T. Let any vertex of T be the root, denoted by v'^. 

We now associate a subset of V{G) with any vertex or edge of T. If v' S V{T) 
is a face, the set s{v') consists of the vertices on the face. If v' is a cut vertex or 
a leaf, s{v') consists of that vertex. If both endpoints of e' G E{T) are vertices 
of G, then s(e') consists of those vertices. Otherwise, let either endpoint of e' is 
a face. Let e' = (fi,W2): ^^^d define s(e') = s{v'i) n s{v2). 

Removing an edge (v' , w') divides T into two components, T^i and T^i , where 
Tyi is the one containing v' . Define B{v' , w') as the subgraph of G corresponding 
to Tyi. Formally, it is the subgraph of G induced by IJu'eyfr ,) s{u'). Thus G is 
divided into two subgraphs B{v',w'), B(w',v'). If v',w' are adjacent faces, the 
subgraphs share two vertices. If v' is a face and w' a cut vertex (or the other way 
around), then the subgraphs share one vertex. Otherwise they do not intersect. 
Define B-{v',w') as the subgraph of G induced by V{G) \ ViB{w',v')). Then 
we can divide G into two nonintersecting subgraphs, B^ [v' ,w') and B{w',v'). 

The following lemma is immediate from the definitions. 

Lemma 12 Given neighbors v' ,w' G V{T), B{v' ,w') and B(w' ,v') are con- 
nected subgraphs of G , and any path from B{v' ,w') to B(w',v') intersects s{{v' ,w')). 
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Given a landmark z, let S{z, W) — {{x, y} £ W x W : d{z, y) ^ d{z, x)} be 
the set of pairs in W^ C V{G) resolved by z. 

Lemma 13 Let e' — {v\w') £ E{T). Assume z e V{B{v' ^w')) and denote 
W — V{B{'w' ,v')). Then S{z^W) ~ S{z,W), where z is the representative of z 
on e' . 

Proof. Let z € V{B{v\ w'j) and x € W. li z = vi or V2, then d{z, x) = d{z, z) + 
d{z,x). If z = We, then d{z,x) = d{z,z) + d{z,x) — 1. So for any x,y £ W, we 
have d{z, x) — d{z, y) = d{z, x) — d{z, y). D 

Lemma 14 Let e' = (t;',w') G E{T) be an edge. If e' corresponds to an edge 
e = (wi,U2) e E{G) with midpoint v^, then S{ve:W) C S{vi,W) U S{v2,W). 

Proof Assume that {x, y} e S{ve, W) \ S{vi, W). So, we have 

d{vi,x) = d{vi,y) and d{ve,x) ^ d{ve,y) 

Without loss of generality, we assume that d{ve,x) < d{Vf,,y). Then 

2 + min{d{vi,x),d{v2,x)) < 5 + inin{d{vi,y),d{v2,y)). 

Putting these equations together implies that d{v2, x) < d{v2, y), that is, {x, y} G 

S{V2,W). ' D 

3.3 Boundary conditions 

Let v' be a dual vertex and p' its parent. We use boundary conditions to describe 
the relation of landmarks in B{v',p') to B{p',v'), and vice versa. Boundary 
conditions can be seen as an abstract data structure which depends on w',p',L 
and satisfies the following: 

1. It consists of two components, one of which depends on landmarks in B{p' , v'), 
and the other on landmarks in B^ (v' ,p'). The components are called upper 
and lower boundary conditions, respectively. 

2. It determines which pairs in B{p' , v') are resolved by landmarks in B~ («', p'), 
and vice versa. 

3. If B{v' ,p') contains landmarks with a representative v on v' , then the bound- 
ary condition specifies which such landmarks occur first and last on a walk 
along the outer boundary of B{v',p') starting at v. 

4. For any v £ s{{p', v')), it specifies whether the set g{v, L) n V{B^{v' ,p')) is 
empty or not. 

5. The number of possible boundary conditions is polynomial. 

The first and the second properties are necessary to be able to run a dynamic 
programming algorithm along T. The third and fourth properties are needed 
when verifying Requirements [2] and [ll respectively. The last property is needed 
to ensure that the algorithm runs in polynomial time. 
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Let us consider the implementation of this data structure. Let e' = (p', v') G 
E{T), s(e') — {vi, U2}. The boundary condition of e' is a tuple 

t = {tu,ti,t„ty^,t,,J e V{{vi,v2,v,})^ X V{V{G)) X {0,1, 2}2. 

The set t„ consists of the representatives oi LnV{B{p' ,v')) on e', and ti consists 
of the representatives of L n V{B~{v',p')) on e'. 

The set tz consists of all the first and last landmarks with given representative 
as property 3 in above requires. Note that i; is technically redundant, since t^ 
determines its value. As an optimization, t^ could be omitted in some cases. 

The value of t^-^ depends on the properties of the set 

Vi = g{vi, LnViB~{v',p'))) n V{B-iv',p')). 

If Vi contains a vertex x which is not a neighbor of V2, then ty-^ = 0. If Vi — \x\ 
with X e M{v\) n A/'(w2), then i„^ = 1. Otherwise V\—% and t^^ — 2. We define 
i^2 similarly, with the roles of v\ and V2 interchanged. 

Lemma 15 The implementation of boundary conditions described above has 
properties 1-5 given above. 

Proof. A boundary condition (t^ ,ti,tz, t^-^ , t^^ ) can be divided into an upper 
boundary condition t„ and a lower boundary condition {ti,tz,ty-^,tv2), so prop- 
erty 1 holds. Property 2 follows from Lemma [TSJ Property 3 follows from the 
definition of t^ . 

To verify that property 4 is satisfied, consider g{vi, L) n V{B^{v' ,p')). This 
set is empty if either t^-^ = 1 and V2 G t^, or if ty-^ — 2. Otherwise, the set is 
nonempty. We can interchange vi and W2, and the same holds. So property 4 is 
satisfied. 

It remains to verify property 5 is satisfied. There are (at most) 0(|y|) edges 
e', and for each such edge there are 2^^ — 1 possible values of t„, 2^^ — 1 possible 
values oiti, 3 possible values of ty^, 3 possible values of ty^, and 0(|F|^) for t^. 
Therefore property 5 holds. D 

3.4 Configurations 

While a boundary condition describes how the landmarks are placed in relation 
to a dual edge e', a configuration describes their relation to a dual vertex v' . 
There are three quite different cases, depending on whether v' is a cut vertex of 
G, a face with two representatives, or a face with at least three representatives. 
The last-mentioned case is quite complicated. 

A configuration L associated with v' is designed to satisfy following: 

1. For any w' e M{v'), it determines which pairs of vertices in B{w' ,v') are 
resolved by landmarks in B{v',w'). 

2. It contains enough information to verify that v' satisfies Requirement ^ 

3. The total number of configurations is polynomial. 
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Essentially, a configuration determines which vertices of s(w') are landmarks, 
and which boundary conditions are allowed for edges {v',w'). 

Implementation of this data structure starts by extending the concept of 
representative to take into account the branch of the dual that contains the 
landmark yielding the representative. Formally, the extended representative of 
z G L on v' is a pair (z,w'). Here z is the representative of z on v' . Moreover, 
if z G s{v'), then w' — v'; otherwise, w' is the neighbor of v' such that z G 
V{B^{w',v')). 

Given a face v' and vertex z € s(w'), there are at most four values of w' such 
that (z, w') is a valid extended representative, namely w', at most two neighboring 
faces, and at most one neighboring cut vertex. For example, in Fig. [91 if z = wig, 
w' can be Fi, F2, F3, or viq. 

We can now define the configurations for each of the three cases. 

Definition 16 Let L be a set of landmarks and v' a dual vertex which corre- 
sponds to a cut vertex v e V{G). Let R consist of all extended representatives of 
L on v' . The configuration associated with v' , L is the pair (w',i?). This is said 
to be a configuration of type I. 

In this case, any extended representative has pattern {v,w') with w' G 
J\f{v') U {w'}. Hence the total number of configurations of type I would be 
2(deg(f )+i)^ However, from here on we only consider the ones that fulfill Re- 
quirement IT] Then there is at most one w' € Af{v') such that {v, w') ^ R. 
Therefore the total number of configurations of type I is 0(|£'(r)|) = 0(|1/|). 

Definition 17 Let L be a set of landmarks and v' a face such that L has two 
representatives on v' . Let R = {zi, z[, Z2, Z2}, where the landmarks are defined 
as in Requirement^ Let the configuration associated with v' and L be the pair 
(«', R). This is said to be a configuration of type II. 

In this case the number of configurations clearly is Odyj'^). 

The remaining case, a face with more than two representatives, is the most 
intricate because of combinatorial explosion: on a face of n vertices, there are 
up to 3" possible sets of representatives. 

Definition 18 Let L be a set of landmarks and v' a face such that L has more 
than two representatives on v' . The configuration associated with v' ,L is a pair 
{v' , R), where R is computed by the (nondeterministic) Algorithm\U This is said 
to be a configuration of type HI. 



As an example of a configuration of type HI, consider Fig. 10 

Note that Algorithm [l] is a part of a mathematical definition. It is not ac- 
tually used in any computations. Algorithm [I] is given a set of landmarks and 
then produces a configuration. In our computations to find a metric base (Al- 
gorithm l4|, we proceed in the opposite direction: we iterate over pairs {v',R), 
and for each of them, we produce a set of landmarks, if possible. The following 
lemma shows that we can restrict to sets R such that \R\ < 20, and furthermore 
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Algorithm 1 Construct-Conf 

Input: Dual vertex «', set L 

1: Let L be the set of representatives of L on v' . 

2: Choose zi,Z2 & L sd maximal mutual distance. 

3: Choose £3 G L \ {zi, 22} which is, if possible, not on a shortest path 21 -^ £2- 

4: Ri -«- {zi,Z2,Z3} 

5: if the shortest path zi -^ £2 is not unique and both paths contain a vertex of 
L \ {zi, Z2} then 

6: Choose 24 G L \ {21, 22} which is on a different shortest path than 23 

7: Ri^RVJ{ii} 

8: end if 

9: for all i,j G {1,2,3} do 
10: if d{zi, Zj) = (s{v') — l)/2 and there is 2„ £ L\R\ which is on the shortest path 

Zi -^ Zj then 
11: i?i ^ i?i U {2„} 

12: end if 
13: end for 

14: Let R consist of all extended representative corresponding to elements of Ri 
15: return R 



the set Q — {z : {z,w') G R for some w'} must satisfy 3 < \Q\ < 5. Because 
\R\ < 20, the total number of configurations of type III that wc need to iterate 
over in our computations thus is OdFp"). 

Lemma 19 Algorithmln produces a set R with \R\ < 20, such that R determines 
the upper boundary condition of (w', v') for any child w' of v' . If v' has parent p' , 
R also determines the component ti of the lower boundary condition of{p',v'). 

Proof. We first show that \R\ < 20. Initially, Ri will consist of three elements. 
Either line [6] adds one more element, or the loop on lines [9HTT] adds at most two 
elements. So \Ri\ < 5, and \R\ < 4|i?i| = 20. 

We now verify the effect of R on the boundary conditions described in the 
lemma statement. Observe that, by Lemma |13[ any two landmarks with the 
same extended representative have the same effect on the boundary conditions 



in Definition 18 Therefore it is indeed justified to talk about the upper/lower 
boundary condition determined by R. 

We first prove the following auxiliary claim. 

Claim 1: Any z G L is on a shortest path Zi ^^ Zj for some i,j G {1, 2,3}. 

Proof: If every representative is in shortest a path zi ^ £2, the claim holds. 
Otherwise, 23 is selected so that it is not in the shortest path zi -^ £2- Then 
shortest path zi -^^ Z2 and zi '^-> £3 intersect only in zi, because otherwise 
d(zi,Z3) > (i(zi,Z2), contradicting the choice of zi and £2- Also shortest paths 
Z2^^ zi, Z2 ^^ £3 intersect only in Z2. Hence the shortest paths cover all vertices 
in s{v'). This proves the claim. D 

We now prove that R determines the upper boundary condition of {w',v') 
for any child w' of v' . So let w' be a child of v' , let e' = {v',w'), and let 
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Fig. 10. Resolving set L = {113,^6, V7,vi5,ui7,wi8} has a configuration {Fi,R) with 
R = {{viSjVie), (via, F2), (fis, F4), {vi,vi)}. Values of the function h are h{v\,F\,L) = 
vi (because g{vi,L) = {vq} C F(B"(«i, Fi))), h{vio,Fi,L) = 0, h{vi:i,Fi,L) = 
h{vie,Fi,L) = F3, h{v2o,Fi,L) = h{v2i,Fi,L) = 0. 



s(e') = {wi,W2}- Let i^ be the upper boundary condition of L on e', and let r^ 
be the upper boundary condition of i? on e'. We need to show that r„ = Z„. 

The inclusion w e r^ holds if and only if there is a (z,w') € -R such that 
u' ^ w;' and z has representative v' on e'. Similarly, v (z tu ii there is z e L with 
such an extended representative. 

Since R corresponds to a subset of L, r„ C t„. To show inclusion in the other 
direction, consider v £ tu- By definition there is z € L with z ^ V{B{w' ^v')) 
that has representative v on e'. Denote the representative of z on u' by z. By 
Claim 1, z is on a shortest path Zi -^ Zj , for some i, j € {1,2,3}. We distinguish 
several cases. 



Case 1: z = 

Since z ^ y(i3(w', w')), z 7^ Ve, and thus w 7^ Ve- If one of Zi,Zj has representative 
V on e' (say Zi), then it must have an extension [zi,u') with m' ^ w' (as in 
Figure 



11 



) Therefore v G r„ and we are done. So suppose that neither has 
representative ij on e'. Then the representative of both Zi and Zj on e' is t^e (the 
midpoint of e). By construction, this implies that {i,j} = {1,2}. In particular, 
one of Zi, Zi equals Ve and the other is antipodal to it (see Fig. 13). Then either 
Z3 or Z4 has representative v on e', and we are done. 

Case 2: z ^ V. 

We distinguish two subcases. 

Case 2a: v — Ve- 
in this case, either one of z^ 
Ve € r„, and we are done. Otherwise, 



equals z, or not. If one of Zi,Zj equals z, then 
= {t;i,U2} and tu = {wi,i'2,We} (see 
Fig. 12 ), and the sets are equivalent in the equivalence relation used for boundary 



conditions. 
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Fig. 11. Illustration for Case 1 in the proof of Lemma [19] 




Fig. 12. Illustration for Case 2a in the proof of Lemma p^ Here v is the midpoint. 



Case 2b: v ^ Ve- 
in this case, Zi and Zj both have representative Ve or u on e'. We again distinguish 
several subcases. 

Case 2b-i: Zi and Zj both have representative v,, on e'. 



In this case, the shortest path Zi 



is not unique. Hence either Z3 or 24 has 



representative v on e' (see Fig. 13), and we are done. 

Case 2b-ii: Both Zi and Zj have representative v on e'. 

In this case, at least one of them (say Zj) has an extension (zj, u') vifith u' 7^ w' 

(see Figure 14 1. Hence w € r„. 



Case 2b-iii: Zi has representative v on e', Zj has representative iig on e'. 
In this case, either Zi = v or not. If Zi = w, then line [TT] of Construct-Conf 
adds an element to Ri — the set L\ Ri is not empty because it contains z — and 
thereby ensures that v g r„. Otherwise, Zi has an extension (zi,u') with u' ^ w' 



(see Figure 15 1. Hence v £ r^. 



This completes the case analysis. It follows that r„ = i„. 

We note that exactly the same arguments hold when w' is not a child of v' 
but its parent p' . D 

Even if both boundary conditions for edges {v' , w[), {v' , W2) are allowed by a 
specific configuration, they may contradict each other. This happens when there 
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Fig. 13. Illustration for Case 1 and Case 2b-i in the proof of Lemma [191 Both Zj and 
Zj have representative Ve on e'. 




Fig. 14. Illustration for Case 2b-ii in the proof of Lemma [TQ) Now Zi 



\s V € s{{v' ,w'i)) n s((w', W2)), and the boundary conditions disagree about the 
value of g{v, L). Hence, the algorithm will only process boundary conditions that 
agree with each other and with the configuration. 

As a tool for avoiding such disagreements, we define a coarser variant of the 
function g(-,-). The function /i : V (G) x V (T) x V {V (T)) -^ V{T)U{9} describes 
which part of the generalized dual contains g{v, L). Let v £ s(v'). Then 

r v' iig{v,L)ns{v')^(!l, 
h{v,v' ,L) — < w' if w' is neighbor of f' and 5(u,L) n V^(i?~(it;',t;')) 7^ 0, 
[ iig{v,L) = %. 

Notice that as long as Requirement [T] holds, the function is well defined. 

Lemma 20 The configuration associated with v' and L determines for any v G 
s{v') whether the equation h{v,v',L) = v' holds. 

Proof. Recall that h{v,v',L) = v' if g{v,L) D s{v') 7^ 0. If the configuration 
(w', R) is of type I, i.e., v' corresponds to a single vertex of G, then h{v, v' , L) = v' 
is not possible. 



Otherwise, u' is a face. By Lemma 13 g{v, L) n s{v') — g{v, L) n s(w'), where 



L is the set of representatives of L onv' . If (u, R) is of type II, then it determines 
L, and we are done. Hence the configuration is of type HI. Let Q = {z : {z, w') G 
R for some w'}. Then Q ^ L. We claim that g{v, L) n s{v') = g{v, Q) n s{v'). 
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Fig. 15. Illustration for Case 2b-iii in the proof of Lemma [19] in case Zi 



Let w,v (z •s(u') be neighbors. If ui ^ g{v,L), then there is z e L with 



d{z, w) < d{z, v). As shown in the proof of Lemma 19 the representative of z on 
v' is on a shortest path zi -^ Z2 for some zi, £2 G Q- Both of them also satisfy 
d{zi,w) < d{zi, v). Hence w ^ g(v, Q). 

li w ^ g{v,Q), there is z G Q with d{z,w) < d{z,v). Replace z with any 
vertex z G L that has representative z on u', and the inequality still holds. 
Hence w ^ g{v,L). D 



3.5 Algorithm 

The algorithm works in a dynamic programming manner, by finding optimal 
resolving sets of subgraphs B{v' ,p') with given boundary conditions. Formally, 
assume that we are given a vertex v' £ V{T) and boundary condition r on the 
edge e' = {v',p'), where p' is the parent of v' . Let m{v' ,t) C V{B~{v' ,p')) be a 
set Lr)V{B~{v' ,p')), where L is a minimal resolving set with boundary condition 
t, if such a set exists. Otherwise m{v',t) — nil. For notational convenience, we 
let |nil| = 00 and nil U A = nil for any A. 

The values of m(v' , t) are computed in a recursive manner: the computation 
of m{v',t) uses the values of m{w',r), where w' is a child of v' . The basic idea 
is to iterate over all configurations on v' . For every configuration, we then find 
an optimal function h and an optimal set of landmarks. 

First, we introduce several subroutines. Algorithm ^ (Intermediate-sol) is 
given a configuration (v',R)^ a boundary condition t, and a function h, and it 
returns an optimal set of landmarks. 

Lemma 21 Algorithm pi runs in polynomial time and correctly computes an 
optimal set of landmarks given a configuration (v'^R), a boundary condition t, 
and a function h. 

Proof. We first consider the complexity of Algorithm [2] The consistency checks 
run in linear time. The for-loop runs over 0(|T^|) neighbors, and for each neigh- 
bor, there are Odl^l"*) potential boundary conditions. Hence it runs in polyno- 
mial time. 
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Next we argue that Algorithm [2] is correct. It returns nil only if a resolving 
set of B{v' ,p') with the given parameters does not exist. The algorithm explicitly 
verifies that Requirement [2] is satisfied. Requirement IT] is verified implicitly by 
constructing a solution which agrees with h. By assumption the Requirements 
are satisfied in subgraphs B{w',v'). Therefore Lemma III] states that a non-NiL 
return value of Algorithm [2] is a resolving set of B{v',p'). Assuming that the 
values of m for all children w' of w' have been computed correctly, the returned 
set will be an optimal set under the conditions given in the lemma statement. D 

Given a configuration (w', R) and a boundary condition t, the next subroutine 



finds an optimal function h. By Lemma 20 the configuration determines whether 



h{v,v',L) = v' holds or not. Also t restricts some values of h. Otherwise, the 
value of h{v, v' , L) may be or w' , where w' is a suitable neighbor of v' , and the 
task is to determine which one of these is optimal. It turns out that the optimum 
can be found by a greedy algorithm (Algorithm [3| . 

Lemma 22 Algorithm[Mruns in polynomial time and returns a smallest resolv- 
ing set of B{v' ,p') that agrees with the parameters. 

Proof. The for- loops run over 0(|V"(r)|) cases. The subroutine Intermediate-sol 
runs in polynomial time. Since |s((w',w'))| < 2, there is a constant number of 
possibilities for k in each iteration. Therefore Opt runs in polynomial time. 

We need to prove that the function h constructed by Opt is optimal. Let h° 
be an optimal function, and define h^ as follows. li v ^ Q after i iterations of the 



loop in fine 10 then h^[v] = h[v]. Otherwise, h^[v] = h°[v]. Therefore ft." = h° 
and /i™ = h, where m = \J^{v')\ - 1. 

First note that changing the value of /i[w] from some other value to increases 
the value of Intcrmediate-sol(w', R, r, h) by at most one. To see this, let L be the 
original return value with g{v, L) — {x}. Then Ln ~ LU {x} is an resolving set 
with g{v,Ln) — 0. 

To conclude the proof, we show that 

|Intermediate-sol(t;', R, r, /i*)| > |Intermediate-sol(w', R, r, h''^^)\. 



Algorithm 2 Intermediate-sol 

Input: Configuration {v',R), boundary condition t, function h 

if the parameters are inconsistent or {v' , R) violates Requirement [2| then 

return nil 
end if 

Initialize W to the set of landmarks on s{v') as described by {v' , R) 
for all w' that are children of v' do 

r <— arg miUj, mlw' , r] such that r agrees with C, h, t 

W <- WUm[w',r] 
end for 
return W nV(B' (v ,p')) 
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Algorithm 3 Opt 

Input: Configuration {v',R), boundary condition t 

Q -^ {Q is the set of vertices for which h[v] is already fixed} 
for all V G s{v') do 

if R or r determine h{v,v',L) then 
h[v] <— the appropriate value 
Q^Qu{v} 

else 

h[v] ^ 

end if 
end for 

for all w' that are children of v' , in clockwise order, starting on the successor of p' 
do 

P^s{(v',w'))\Q 

Find k that minimizes |Intermediate-sol(u', i?, r, fc)| such that 

h[v] and k[v] differ only for v £ P. li possible, choose a solution with k[vi] — for 

the last (in clockwise order) element vi £ P 

h -i^ k 

Q^QU{v: h[v] /0} 
end for 
return Intermediate-sol(i;', R, r, h) 



Let vi be as in the algorithm and consider the {i + l)-th iteration of the second 
for-loop in Opt. Observe that h^ and /i*"*"^ differ only in s{{w',v')) \ Q. If they 
differ in V2 7^ vi, then h'^[v2] = and changing the value to /i*"''^[v2] does not 
increase the target value. Then, if h'^{vi) 7^ h'^^-^{vi), we first replace h^{vi) with 
and then with h^'^^{vi). By the argument above, the first change may increase 
the size by at most one. The second change must decrease the size, because 



otherwise the optimization made in line 12 would not change h[vi]. Therefore 



the changes do not increase the size of the corresponding resolving set. D 

Finally, we are ready to present the main algorithm (Algorithm [4]) and its 
correctness proof. 

Theorem 23 AlgorithwY^ correctly computes the values ofm[v',Y\, and returns 
a metric base in polynomial time. 



Proof. The algorithm runs in polynomial time, because there is a polynomial 
number of dual vertices, boundary conditions, and configurations. 

We prove by induction that m[?;',r] is computed correctly. Assume that the 
values of m[w',t] have been computed correctly for any child w' of v' . Then, 
because Opt (Algorithmlsl) works correctly, the value oim[v' , r] will be computed 
correctly. Similarly, the return value will be computed correctly, since Opt works 
correctly and uses correct values of m[w', r]. D 
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Algorithm 4 Metric-D 



Input: Graph G, its generahzed dual tree T 

for all v' G V{T) \ {v'^}, boundary condition r do {recall that v^ is the root} 

m[v',r] <— NIL 
end for 

for all v' eV{T)\ {v'r}, in postorder do 
p' -f- the parent of v' 
for all configuration {v',R) do 

for all boundary condition condition r on (v' ,p') do 
if |Opt(v',i?,r)| < |m[u',r]| then 

m[v\r] 4- Opt{v\R,r) 
end if 
end for 
end for 
end for 
W ^ V{G) 

for all configuration {v'^,R) do 
if |Opt(t;;,7?,NiL)| < \W\ then 

W -^ Opt«,i?,NIL) 

end if 
end for 
return W 



4 Conclusions and Open Problems 

We have proved that Metric Dimension is NP-hard for planar graphs, even 
when the graph has bounded degree. We also gave a polynomial-time algorithm 
to solve the problem on outerplanar graphs. Our algorithm is based on an inno- 
vative use of dynamic programming which allows us to deal with global problems 
(in the sense described in the introduction), and which may be amendable to 
other such problems. One possible candidate is the identifying codes problem. 
Let Nr{v) = {u e V{G) \ d{u,v) < r}. Then given as input a graph G and an 
integer r > 0, the identifying codes problem is to find a smallest D C V{G) 
such that Nr{v)nD y^ H) for all v e V{G), and Nr{v) n D ^ Nriu) n D for 
all distinct u,v E V{G). This problem is known to be NP-hard for general and 
planar graphs, and to have polynomial-time algorithms for cycles, paths, and 
trees (see Jj and references therein). We believe that the techniques developed 
in the present paper can yield a polynomial-time algorithm for the identifying 
codes problem on outerplanar graphs when r is part of the input. 

We also pose some open problems for Metric Dimension. First, it would be 
nice to extend our results to fc-outerplanar graphs. Unfortunately, the technique 
described in the present paper does not seem to generalize in this direction, 
as edges and cut vertices are no longer convex separators in this case. Even if 
the problem turns out to be solvable on fc-outerplanar graphs by a polynomial- 
time algorithm, it is not clear that such an algorithm could be used to derive a 
polynomial-time approximation scheme for Planar Metric Dimension. The 
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quest for such an approximation scheme or even for a constant approximation 
algorithm is an interesting chaUenge in its own right. 

Another interesting Une of research is the parameterized complexity of Met- 
ric Dimension. Daniel Lokshtanov [17^ posed this problem at a recent Dagstuhl 
seminar on parametrized complexity. Moreover, he conjectured that the problem 
could be VF[l]-complete. We hope that our results can serve as a starting point 
for further study in this direction. 

Acknowledgment: The authors thank David Johnson for sending a copy of the 
NP-completeness proof. 
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